Magnetic properties of the S=l/2 quasi square lattice antiferromagnet 
CuF 2 (H 2 0)2(pyz) (pyz=pyrazine) investigated by neutron scattering 
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We have performed elastic and inelastic neutron scattering experiments on single crystal samples 
of the coordination polymer compound CuF2(H20)2(pyz) (pyz=pyrazine) to study the magnetic 
structure and excitations. The elastic neutron diffraction measurements indicate a collinear antifer- 
romagnetic structure with moments oriented along the [0.7 1] real-space direction and an ordered 
moment of 0.60 ± 0.03 /xs/Cu. This value is significantly smaller than the single ion magnetic 
moment, reflecting the presence of strong quantum fluctuations. The spin wave dispersion from 
magnetic zone center to the zone boundary points (0.5 1.5 0) and (0.5 1.5) can be described by a 
two dimensional Heisenberg model with a nearest neighbor magnetic exchange constant J2d = 0.934 
± 0.0025 meV. The inter-layer interaction J perp in this compound is less than 1.5% of J^d- The spin 
excitation energy at the (0.5 0.5 0.5) zone boundary point is reduced when compared to the (0.5 
1 0.5) zone boundary point by ~ 10.3 ± 1.4 %. This zone boundary dispersion is consistent with 
quantum Monte Carlo and series expansion calculations for the S=l/2 Heisenberg square lattice 
antiferromagnet which include corrections for quantum fluctuations to linear spin wave theory. 

PACS numbers: 75.30.Ds, 75.40Gb, 75.50.Ee 



I. INTRODUCTION 

Metal-organic systems with a 3d 9 electron configura- 
tion, such as in Cu 2+ , are expected to undergo a Jahn- 
Teller distortion. In octahedral coordination, this typi- 
cally elongates one axis of the octahedron and removes 
the degeneracy of the e g orbitals d x 2_ y 2 and d 2 2. This 
effect widely occurs in molecular systems as well as in 
other materials such as the colossal magnetoresistancc 
manganites. 1,2 The Jahn- Teller distortion is extremely 
sensitive to bond distances and as such applied pressure 
can strongly influence or even induce a Jahn- Teller dis- 
tortion. Depending on the elongation axis in the crystal 
structure, the pattern of orbital overlaps and related ex- 
change interactions can vary significantly in turn leading 
to disparate magnetic ground states. Finding materials 
where a Jahn- Teller distortion can be tuned to act as a 
magnetic switch has potential for applications in techno- 
logical devices and, thus, is of great interest. 

The copper-based coordination polymer magnet 
CuF2(H20)2(pyz) (pyz=pyrazine) appears to be a 
model material for studying the switching of magnetic 
properties due to changes in the Jahn- Teller axis. 3,4 
CuF2(H20)2(pyz) crystallizes in a monoclinic structure 
(space group P2i/c) with a=7.6926 A , b=7.5568A , 
c=6.897A and /3=111.065° under ambient conditions. 5,6 
The structure consists of CUF2O2N2 octahedra which 
form a 2d network in the be plane. The magnetic sus- 
ceptibility shows a broad peak at 10.5 K consistent with 
two dimensional (2d) short range correlations and three 
dimensional long-range order below Tn = 2.6 K due to 
weak coupling along the a-axis. At ambient pressure, 
the Jahn- Teller axis is along the N-Cu-N axis (a direc- 



tion). Hence the A x i_ y i magnetic orbital lies in the be 
plane to form a two dimensional (2d) antiferromagnetic 
(AFM) quasi square lattice where the near neighbor dis- 
tances are the same, but the interior angles deviate by 
± ^5°. First principles electronic structure calculations 
confirm this configuration. 5 These calculations give an 
intra-plane AFM exchange interaction of about 13 to 19 
K and an inter-plane AFM exchange interaction of only 
1% of the intra-plane exchange interaction. An impor- 
tant feature of CuF2(H20)2(pyz) is that under an ap- 
plied pressure of 0.9 GPa, a switch of the Jahn- Teller 
axis from the N-Cu-N to the O-Cu-0 (c direction) bond 
occurs. 3 When pressure increases to 3.1 Gpa, the Jahn- 
Teller axis switches again from O-Cu-0 bond to F-Cu-F 
bond (b direction). Correspondingly the magnetic inter- 
actions are expected to vary at different pressures: from 
the ambient pressure 2d quasi square lattice (exchange 
path Cu-F- • H-O-Cu) to one dimensional (Id) chain in- 
teractions (exchange path Cu-pyz-Cu) at 0.9 GPa. In- 
deed, the magnetic susceptibility at ambient pressure can 
be well described by the 2 d Heisenberg square lattice 5 
while above 0.9 GPa the magnetic susceptibility exhibits 
a broad peak consistent with Id magnetic correlations 3 . 
Therefore a microscopic understanding of the magnetic 
properties at ambient pressure would be an illuminating 
step towards a complete understanding of the interesting 
pressure dependent behavior of CuF 2 (H 2 0)2(pyz). 

In this paper we present the results of a study 
of deuterated CuF 2 (H 2 0) 2 (pyz) (CuF 2 (D20) 2 (d4-pyz)) 
with elastic and inelastic neutron scattering under ambi- 
ent pressure. The neutron diffraction results show that 
the sample has a collinear AFM structure with the mo- 
ments lying in ac-plane along the [0.7 1] direction. The 



2 



spin wave dispersion extracted from the inelastic neutron 
spectra can be well described by a nearest-neighbor 2d 
Heisenberg model. The spin wave dispersion found along 
the inter-layer direction could not be observed within the 
instrumental resolution of FWHM=0.034 meV. This in- 
dicates that the inter-layer exchange interaction J per p is 
weak as expected for a 2d system where the spin wave 
dispersion does not depend on the out of plane direction. 
Consequently, CuF 2 (H 2 0)2(pyz) is a good example of 
quasi-2d system where the spin wave dispersion depends 
exclusively on K and L. Within the 2d magnetic plane, 
the dispersion along the zone boundary points has been 
measured and is in accord with prediction of the quan- 
tum Monte Carlo and series expansion calculations 13 ' 15,17 
for quantum fluctuation corrections to linear spin wave 
theory for the S=l/2 square lattice Heisenberg antiferro- 
magnet. 



II. EXPERIMENTAL DETAILS 

In order to minimize attenuation of the beam caused 
by the large incoherent scattering cross-section of hydro- 
gen, fully deuterated CuF2(D20)2(d4-pyz) crystals were 
synthesized for this study. Substitution of deuterium for 
hydrogen does not alter the crystal symmetry and in- 
duces only small changes in the lattice parameters. 6 To 
ensure that structural parameters of the samples stud- 
ied here are in accord with previously published values 6 
neutron diffraction data (not shown) were collected on 
the four circle neutron diffractometer HB-3A at the High 
Flux Isotope Reactor (HFIR), Oak Ridge National Lab- 
oratory. The magnetic structure was studied with the 
thermal triple- axis spectrometer HB-1 at the HFIR. A 
crystal with a mass of 0.01 g was studied using HB-3A 
while a larger crystal with mass 0.1 g was studied with 
HB-1. A silicon monochromator and no analyzer were 
used for the measurements on HB-3A and unless other- 
wise noted all triple-axis measurements were performed 
using pyrolytic graphite (002) monochromator and ana- 
lyzer crystals. 

Inelastic neutron scattering experiments were per- 
formed using the cold neutron triple-axis spectrometer 
CG-4C (HFIR), the thermal triple-axis HB-1 A (HFIR) 
and the cold neutron chopper spectrometer CNCS 7 at 
the Spallation Neutron Source (SNS), Oak Ridge Na- 
tional Laboratory. The inelastic neutron scattering ex- 
periments were performed on four co-aligned crystals 
with total mass of 0.85 g and total mosaic of 0.7°. For 
the CG-4C experiment, measurements were made with 
fixed final energy, Ef , of 5 meV and 3.7 meV resulting in 
an elastic energy resolution of about 0.3 meV and 0.17 
meV, respectively. For the HB-1A experiment, the inci- 
dent energy Ei was fixed at 14.7 meV with a pyrolytic 
graphite (002) monochromator and a beryllium (002) an- 
alyzer was used to obtain an elastic energy resolution of 
about 0.58 meV. For the CNCS measurements, the inci- 
dent energy Ei was fixed at 3 meV and 1.5 meV to allow 
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FIG. 1. (a) Temperature dependence of magnetic Bragg 
peaks (a): (0.5 1 0) and (b): (-0.5 1) for CuF 2 (D 2 0)2(d4- 
pyz). Insets display the temperature dependent peak inten- 
sity. Solid lines are guides to the eye. The data were collected 
on the HB-1 triple-axis instrument in HK0 and H0L scatter- 
ing planes, respectively. 

examination of the entire spin wave spectrum. The en- 
ergy resolution (FWHM) was about 0.08 meV and 0.034 
meV at the elastic position for Ei — 3 meV and 1.5 meV, 
respectively. 



III. RESULTS AND DISCUSSION 

A. magnetic structure 

Temperature dependent neutron diffraction reveal the 
presence of additional peaks for temperatures less than 
2.5 K at wave vectors Q=(0.5 K L) with integer values 
of K and L and where K+L is an odd number. This in- 
dicates that in real space the magnetic lattice is doubled 
along the a direction compared to nuclear lattice while 
in the bc-plane the nuclear lattice and magnetic lattice 
remain the same size. Fig. 1 shows the temperature 
dependent magnetic Bragg peaks at both (0.5 1 0) and 
(-0.5 1) collected on the HB-1 instrument using the fi- 
nal energy £y=13.5 meV(A=2.4617 A). In the insets the 
temperature dependent peak intensity for (0.5 1 0) and 
(-0.5 1) are presented. A power-law fit of the form (1- 
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gives a Tn of 2.51 ± 0.06 K and a value of (3 
of 0.24 ± 0.07. This low value of /3 is likely consistent 
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calculated intensity 

FIG. 2. The observed intensity plotted as a function of 
calculated intensity for CuF2(D20)2(d4-pyz) at T=0.25 K. 
Nuclear and magnetic reflections are denoted by diamonds 
and circles, respectively. The solid line is a guide to the eye. 
The data were collected on HB-f in different scattering planes. 
The insets depict different views of the unit cell. The arrows 
on Cu atoms indicates the magnetic moment configuration. 
The 2d magnetic behavior of CuF2(D20)2(d4-pyz) originates 
in the be plane. 

with 2d critical behavior, however, care should be used 
when interpreting this finding as a more rigorous series 
of measurements is required to definitively establish the 
critical exponents. The value of Tjv determined from the 
powerlaw fit is slightly smaller than the value reported by 
muon spin relaxation measurements of a non-deuterated 
sample 5,8 but is in reasonable agreement with the value 
of Tjv= 2.53 K determined from magnetic susceptibility 
measurements of fully deuterated samples. 9 

A total of 26 nuclear and 30 magnetic reflections 
have been collected in HKO and HOL scattering planes 
at T=0.3 K using HB-1. Representation analysis for 
the general case of propagation vector (0.5 0) in the 
structure symmetry P21/c have been performed using 
SARAh- Representational Analysis. 10 The results show 
that there are two possible irreducible representations 
as shown in table I. Using the basis vectors in table I, 
assuming the magnetic moment direction M=(M X M y 
M z ) varies freely, the magnetic reflections can be fitted. 
For both cases, when M y ^ 0, the spins will configure 
in noncollinear way. Fitting of moment direction was 
carried out by minimizing 

v 2 = — V[ — " — V (^)l 2 (1) 
A AT / ; L pobs AT / ^ pobs II \ > 

where N is the number of reflections. Fits considering 
both representations were performed and the irreducible 
representation denoted by Ti yielded the best agreement 
with the data. The resulting collinear spin structure has 
moments along the [0.7 1] real space direction. For sim- 
plicity, the dipole approximation to the Cu 2+ form factor 
has been used. This may lead to a systematic error in 



both the direction and magnitude of the magnetic mo- 
ment due to the anisotropy of the d x 2_ y 2 orbital. 11 We 
note that the Ti irreducible representation with 1^1^=0 
results in an antiparallel alignment of the spins on Cu 
site 1 (atomic coordinate (0 0)) and Cu site 2 (atomic 
coordinate (0 0.5 0.5)) consistent with the observed prop- 
agation vector which requires K+L=odd. The resulting 
magnetic structure for CuF2(H20)2(pyz), shown in Fig. 
2, is G-type (nearest neighbor) antiferromagnetic order. 

TABLE I. Basis vectors for Cu site 1 (atomic coordinates (0 
0)) and Cu site 2 (atomic coordinates (0 0.5 0.5)) determined 
from the representational analysis 10 for space group No. 14 
(P2i/c) and magnetic propagation vector (0.5 0). Here 
IR represents irreducible representations, BV represents basis 
vectors. 
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Interestingly the moment direction [0.7 1] in real 
space does not appear to be along an obvious structural 
direction. The selection of a specific moment direction 
suggests some anisotropy in the spin Hamiltonian. 

The absolute value of the magnetic moment M was de- 
termined by comparing the intensity of magnetic and nu- 
clear Bragg reflections. To carefully estimate the moment 
size, the instrument resolution was taken into account 
using Reslib 12 and the Debye- Waller factor has been in- 
cluded. Our calculations yield an ordered moment of 0.60 
± 0.03 hb/Cu. This value is much smaller than the spin- 
i free ion moment. In low dimensional systems a mo- 
ment reduction occurs due to quantum fluctuations and 
has been observed in a number of materials. 22,23 The ob- 
served and calculated intensities of magnetic and nuclear 
reflections together with the resulting spin arrangement 
are presented in Fig. 2. 

B. spin dynamics 

The spin wave dispersion has been extracted from in- 
elastic neutron scattering measurements. We first dis- 
cuss inelastic neutron scattering data collected using the 
triple-axis spectrometer CG-4C in the HKO scattering 
plane. Both constant-Q and constant-E 1 scans in the 
HKO scattering plane have been performed. Constant- 
E scans at several different energy transfer are plotted in 
Fig. 3(a) and constant-Q scans are shown in Fig. 3(b). 
These measurements were carried out along the (0.5 K 0) 
direction where K=l corresponds to the magnetic zone 
center. The excitation energy at the zone boundary point 
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FIG. 3. (a) Constant-i? scans along (0.5 K 0) collected 
in the HKO scattering plane at T = 1.5 K by the CG-4C 
spectrometer, (b) Constant-Q scans collected in the HKO 
scattering plane at T = 1.5 K by the CG-4C spectrometer. 
The solid lines in both panels are fits with Gaussians. 

Q=(0.5 1.5 0) is E zbK =2.24±0.04 meV. Within instru- 
mental resolution (FWHM=0.17±0.01 meV), no energy 
gap at the antiferromagnetic zone center has been ob- 
served. 

In Figs. 4(a), (b) and (c) the spin wave dynam- 
ics have been studied with the sample aligned in the 
H0L scattering plane with the time-of-flight spectrome- 
ter CNCS. Fig. 4(a) shows that the spin wave excitation 
is essentially independent of H at the excitation energy 
Ai5=0.25 meV, implying only weak interactions along 
H . In Fig. 4(b), the full dispersion along the L-direction 
is shown. The zone boundary energy in this direction 
is i?2bi=2.204±0.003 meV. This value is consistent with 
the previously determined zone boundary energy along 
the indirection (E z i,k)- To search for a small energy 
gap, measurements were performed on CNCS using an 
incident energy of 1.5 meV resulting in an energy resolu- 
tion of FWHM=0.034±0.0007 meV). The measurements 
were centered on the (0.5 1) magnetic zone center. The 
data obtained along the L-direction is shown in Fig. 4(c). 
Again, within instrumental resolution (FWHM=0.034 
meV), no energy gap in L dispersion curve was observed. 
Furthermore, no dispersion was observed along the H- 
direction for L=l. We note that previous DFT calcu- 
lations predict a inter-layer exchange constant of 0.011 
meV 5 and the resulting interlayer dispersion could not be 
observed with our instrumental resolution. Nearly iden- 
tical exchange interactions along the K and L directions 
together with the lack of dispersion along the ii-direction 
indicates the magnetism in CuF2(D20)2(d4-pyz) is pre- 
dominately 2d. 

The inter-layer TJ-direction magnetic exchange was 
also probed with time-of-flight spectrometer CNCS in a 
high resolution mode with Ei — 1.5 meV. The H vs. E 



slice, for L=l (not shown here) is featureless and no obvi- 
ous spin excitation has been observed within instrument 
resolution. Due to the fact that in spin wave theory the 
zone boundary energy of the spin wave excitation is of 
order 2 J, our data show that the inter-layer exchange in- 
teraction Jp er p is expected to be smaller than 0, ° 34 meV 
= 0.017 meV. Considering the i^ fcL =2.204 meV which 
gives the intra-layer J of 1.102 meV, it is reasonable to 
estimate that the inter-layer exchange interaction is less 
than 1.5% of the intra-layer exchange interaction. This 
ratio is in agreement with that determined by both DFT 
calculations and by fitting a Heisenberg square lattice 
model to the magnetic susceptibility which in both cases 
yields J pe r P /J2d ~ 1%- 5 

Since no energy gap has been observed in the 2d square 
lattice plane at the magnetic zone center within instru- 
mental resolution it is reasonable to believe that any 
spin exchange anisotropy must be very small. How- 
ever, the magnetic moments select the [0.7 1] direc- 
tion suggesting nevertheless that a spin anisotropy is 
present. As Cu 2+ moments should have no single ion 
anisotropy and the presence of inversion symmetry pro- 
hibits Dzyaloshinskii-Moriya interactions, we have in- 
cluded exchange anisotropy of the following form: 

H = [SI • S; + A(S? • S* + Sf ■ Sf)} (2) 

where the summation is over the nearest neighbors in the 
2d plane, J is the effective intra-layer exchange parameter 
and A is the exchange anisotropy parameter. 

Linear spin- wave theory yields the spin wave dispersion 

hw Q = 2Jy/l - A 2 cos 2 (Ktt) cos 2 (Ltt) (3) 

where Q is a function of K and L. The dashed lines in the 
inset of Fig. 6 indicate the diamond pattern on which the 
Cu-ions are positioned. Under the approximation that 
the diamond is replaced by a square, the abscissa labels 
in Fig. 6 are related to square lattice notation in the fol- 
lowing way: (0.5 1 0) and (0.5 1) correspond to (tt tt); 
(0.5 1.5 0) and (0.5 1.5) to (tt/2 tt/2)); and finally the 
point (0 0.5 -0.5) to (tt 0). In Fig. 6 we plot the extracted 
dispersion relation for CuF2(D20)2(d4-pyz). The exper- 
imental data in the 2d plane were determined from com- 
binations of triple-axis and time-of-flight measurements. 
The solid line in Fig. 6 is the expected dispersion behav- 
ior of the classical Heisenberg model (A = 1) with J= 
1.102 ± 0.003 meV. 

In the S=i Heisenberg square lattice antiferromagnet, 
it is well established that the spin dynamics can be well 
described by the classical linear spin wave theory with the 
inclusion of quantum corrections. 13-18 The net effect of 
quantum corrections is an overall renormalization factor, 
Z c ss 1.18, resulting in an effective coupling constant J e // 
= Z c J2d when considering the dispersion from the (tt tt) 
zone center to the (7r/2 7r/2) zone boundary. Thus, our 
extracted exchange constant, J, is an effective coupling 
constant and the resulting J 2( j is 0.934 ± 0.0025 meV 
within the square lattice approximation. This value is 
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FIG. 4. Inelastic neutron scattering data from the time-of-flight spectrometer CNCS. (a) Contour map of intensity in the 
HOL scattering plane with A_E=0.25 meV and T = 1.5 K. Background scattering determined from an empty sample holder 
measurement has been subtracted, (b) The spin wave dispersion along the (1.5 L) direction. The data were collected at T = 
1.5 K with Ei = 3 meV. Background scattering determined from an empty sample holder measurement has been subtracted, 
(c) Low energy spin wave dispersion along the (0 L) direction in high resolution mode at T = 1.5 K with Ei = 1.5 meV. 
Intensities in all panels are given in arbitrary units. The integration range is indicated along with the specified direction at the 
top of each panel. Note that the dispersion was found to be independent of the value of H an as such the integration range 
over H was chosen to be large to increase statistics. 
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FIG. 5. Constant-Q scans along the zone boundary points 
ranging from (0 1 -0.5) to (0 0.5 -0.5). The solid line are fits 
to a Gaussian, the data were collected on HB-1A. 

consistent with estimates from the magnetic susceptibil- 
ity and magnetization measurements of fully deuterated 
samples where Jid is about 0.94 meV. 9 

The dispersion resulting from the anisotropic exchange 
yields a zone center energy gap of E zc = 2JV1 — A 2 . As 
noted previously, we have not observed such a gap within 
the instrumental resolution of 0.034 meV. This allows us 
to place a lower bound on A. The coupling constant 
extracted from the zone boundary measurements is J— 
1.102 ± 0.003 meV. Placing an upper bound on the gap 
energy of 0.034 meV results in a value of A (the ratio 
Jxy/Jz) of at least 0.99988. The result would be a very 
small anisotropy resulting in a slightly larger J z which 



would favor spin orientation along the z-axis. Note that 
the z-axis is the spin wave quantization direction, in this 
case the real-space [0.7 1] direction. This same spin 
anisotropy may be related to the observed low-field spin- 
flop like transition observed for CuF2(H2 0)2(pyz) 5 

The spin wave dispersion has been examined along the 
antiferromagnetic zone boundary. In Fig. 5 we present 
constant-Q scans along (0 K -0.5) from K=l ( square 
lattice (ir/2 ir/2)) to K=0.5 (square lattice (ir 0)). In 
Fig. 6 the observed dispersion between (0 1 -0.5) and (0 
0.5 -0.5) has been plotted. Compared to (0 1 -0.5), the 
excitation energy at (0 0.5 -0.5) has been suppressed by 
10.3±1.4%. This observation is inconsistent with linear 
spin wave theory where there is no dispersion along the 
zone boundary. However, proper inclusion of quantum 
fluctuations using series expansion and quantum Monte 
Carlo techniques for the S=l/2 2d AFM square lattice 
predict a dispersion along the zone boundary with the en- 
ergy at (tt 0) 7-10% lower than that at (tt/2 tt/2). 13 < 15 ' 17 
Experimental behavior consistent with these calculations 
have been seen in several model systems. 19-22 However, 
there are some materials where conflicting experimental 
observations exist. For instance, in La2Cu04 21 the en- 
ergy scale at (tt 0) is about 7% larger than the value at 
(7r/2 7t/2) in contrast to the expected quantum behavior 
which has been attributed to ring exchange among four 
spins in the Cu02 plane. In the coordination polymer 
compound Cu(pyz)2 (004)2, 22 the excitation energy at 
(tt 0) is reduced by 11.5% compared to the energy at 
(tt/2 tt/2). The authors believe this value is larger than 
the expectation of series expansion and quantum Monte 
Carlo calculations and they attribute this stronger sup- 
pression to the next-nearest-neighbor interactions which 
enhances quantum fluctuations. In K2V3O8 23 a strik- 
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FIG. 6. Summary of the spin wave dispersion for 

CuF 2 (D 2 0)2(d4-pyz). The data were collected using CNCS 
(solid circles), CG-4C (solid squares) and HB-1A (solid trian- 
gles). The solid line represents the dispersion of Heisenberg 
linear spin wave theory with a nearest neighbor interaction. 
The dash-dotted line represents the results of series expan- 
sion to higher order. 15 The inset depicts the bc-plane. The 
dashed lines show the square lattice with a nearest neighbor 
exchange interaction. 

ing feature of two modes near the zone boundary point 
(tt/2 it/ 2) has been observed experimentally and no fully 
satisfactory explanation has yet been given for this be- 
havior. For our measurements on CuF2(D2 0)2(d4-pyz), 
in Fig. 6 the dash-dotted line is the results of calcula- 
tions including quantum fluctuations with nearest neigh- 
bor interactions 15 . It is clear that the data points follow 
the expected theory and, therefore, we interpret the zone 
boundary dispersion behavior in CuF2(D20)2(d4-pyz) as 
being consistent with the S=l/2 near neighbor 2d AFM 
square lattice with inclusion of quantum corrections. 

IV. CONCLUSION 

We have performed a series of elastic and inelastic 
neutron scattering experiments to study the 2d spin- 



\ quasi square lattice antiferromagnet CuF2(D20)2(d4- 
pyz). The 2d magnetic lattice in the bc-plane is the same 
size as nuclear structure lattice while along the a direc- 
tion the magnetic lattice is doubled. The ordered mag- 
netic moment of the Cu 2+ ions is found to be 0.60 ± 
0.03 /is which is significantly reduced from the expected 
1 /is, indicating strong quantum fluctuations. The spins 
adopt a collinear antiferromagnetic alignment for nearest 
neighbor sites in the 2d plane and are oriented along the 
real-space [0.7 1] direction. 

The spin dynamics of CuF2(D20)2(d4-pyz) have been 
studied using inelastic neutron scattering data obtained 
from both triple-axis and time-of-flight instruments. The 
2d spin wave dispersion can be described by a nearly 
isotropic spin Hamiltonian with small nearest neighbor 
interaction Jid = 0.934 ± 0.0025 meV and very weak 
inter-layer coupling. Along the magnetic zone boundary, 
the excitation energy shows a 10.3±1.4% dispersion con- 
sistent with square lattice calculations including quan- 
tum fluctuations. 
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